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ABSTRACT 


A  method  which  is  based  on  Lagally's  theorem  is  presented  for  computing 
the  resistance  augmentation,  or  thrust  deduction,  for  submerged  hydrofoil-propeller 
arrangements.  The  effect  on  the  magnitude  of  thrust  deduction  of  foil  lift,  propel¬ 
ler  load,  foil  span,  foil  flap  angle,  and  propeller-foil  spacing  was  investigated. 

Computed  examples  and  experimental  thrust-deduction  results  for  several  hypo¬ 
thetical  arrangements  showed  good  agreement. 

INTRODUCTION 

In  recent  years,  a  great  deal  of  theoretical  and  experimental  work  has  been  performed 
in  the  field  of  propulsion  interaction.  4“ 8  In  particular,  considerable  emphasis  has  been  placed 
on  propulsion  interaction  effects  as  associated  with  submerged  bodies  of  revolution  with  ap¬ 
plication  to  submarine  powering.  As  a  result,  important  data  and  criteria  are  now  available  for 
determining  the  coefficients  of  interaction,  i.e.,  wake  fraction  and  thrust  deduction,  for  various 
hull  and  propeller  arrangements.  In  contrast,  this  report  deals  with  propulsion  interaction  phe¬ 
nomena  for  hydrofoil  and  propeller  combinations  with  application  to  powering  of  high-speed, 
hydrofoil-supported  craft.  Although  the  investigations  already  mentioned  have  used  singularity 
systems  which  are  similar  to  those  used  herein,  an  application  leading  to  design  criteria  for 
hydrofoils  has  not  been  done  previously.  The  present  work  differs  from  solutions  for  a  body  of 
revolution  principally  in  the  detail  consideration  of  the  singularities  which  represent  the  pro¬ 
peller.  The  propeller  singularities  we  based  on  total  flue fu  Air*:  thrust  as  obtained  from  quasi¬ 
steady  propeller  theory  using  empirical  viscous-wake  dahs-,. 

The  study  reported  herein  proposes  to  establish  an  analytical  irethc*!  for  estimating  the 
resistance  augmentation  due  to  propeller  action  on  submerged  norea  /ifcating  hydrofoil  configu¬ 
rations.  An  application  of  Lagally’s  theorem  is  used  to  find  the  uu-.jt.  f.fced  resistance  of  a 
hydrofoil;  this  application  gives  solutions  for  (1)  the  effect  of  various  parameters  on  the  aug¬ 
mented  resistance  of  hydrofoils  and  (2)  propeller  performance  in  a  nonuniform  wake. 

To  verily  the  analytical  study,  experiments  were  made  on  various  propeller-hydrofoil 
arrangements  in  the  24-in.  water  tunnel  at  the  David  Taylor  Model  Basin.  The  scope  of  the 
experimental  work  included: 

1.  Determination  of  the  augmented  resistance  of  a  lifting  foil  for  conditions  which  repre¬ 
sent,  with  and  without  nacelle,  infinite  span  and  semi-infinite  span. 

2.  Investigation  of  the  effect  of  variations  in  propeller  load,  propeller  spacing,  and  foil 
angle-of-attack  on  the  augmented  resistance  of  the  foil. 

In  this  report,  the  theoretical  equations  which  provide  a  solution  of  the  interaction 
force  are  discussed  first,  and  then  necessary  expressions  and  data  for  the  wake  of  the  foil 

^Rihmcn  m  listed  on  page  36. 


and  total  propeller  thrust  are  derived.  Finally,  the  computational  results  and  the  experimental 
phases  of  the  work  are  presented  and  compared. 

This  work  was  accomplished  with  the  support  of  the  Bureau  of  Ships  Fundamental  . 
Hydromechanics  Research  Program. 

STATEMENT  OF  PROBLEM 

In  potential  flow,  mathematical  singularities  which  represent  solid  bodies  or  boundaries 
are  widely  used.  The  Lagally  steady-motion  formula5,7  relates  the  force  experienced  by  a 
source  of  given  output  (or  a  sink  of  given  input)  in  an  arbitrarily  formed  flout  to  the  velocity  . 
which  the  surrounding  flow,  undisturbed  by  the  source,  possesses  st  the  location  of  die  source. 
Lagally’s  steady-motion  theorem  can  be  used  to  obtain  a  solution  of  thei  quasi-steady  inter¬ 
action  force  (thrust  deduction)  for  a  propeller-hydrofoil  system.  The  thrust-deduction  problem 
is  analogous  to  the  familiar  problem  of  finding  the  forces  between  two  bodies  in  n  flow  field. 

A  force  arises  from  the  mutual  influence  of  the  foil  on  the  propeller  and  the  propeller  on  the 
foil.  Simply  stated,  a  foil  affects  propeller  performance  through  its  wake,  and  a  foil  (from  the 
point  of  view  of  resistance)  is  affected  by  the  induced  velocity  field  of  a  propeller.  .  . 

Lagally’s  theorem  provides  a  means  of  circumventing  the  detail  of  integrating  pressures 
over  the  foil  surface  to  obtain  the  thrust-deduction  force.  Since  the  interaction  force  is  to  be 
obtained  directly  from  relations  between  singularities,  the  main  problem  is  to  find  appropriate 
generating  singularities  for  a  hydrofoil  and  a  propeller.  Essential  parts  of  ths  problem  are 
determination  of  the  disturbing  wake  velocities  due  to  a  foil  situated  in  a  uniform  potential 
flow  and  determination  of  propeller-fluctuating  thrust*  Specifically,  these  essential  parts  entsr 
into  the  solution  of  the  interaction  or  thrust-deduction  force,  using  the  Lagally  theorem,  by 
providing  the  required  foil-disturbance  velocity  and  propeller-sink  input.  Propeller-thrust 
fluctuations  were  estimated  by  calculating  quasi-steady  propeller  forces.  Inasmuch  as  pro¬ 
peller-fluctuating  thrust  is  caused  mostly  by  the  nonuniformities  (as  seen  by  the  propeller) 
in  the  wake  of  the  hydrofoil  and  the  viscous  part  is  large,  expressions  for  the  viscous  wake 
must  also  be  obtained.  For  this  purpose,  NACA  empirical  wake  data  were  used  ss. a  basis 
for  obtaining  viscous-wake  distributions. 

The  principal  assumptions  snd  limitations  that  are  involved  in  the  analysis  and  treat¬ 
ment  of  the  problem  are  summarized  as  follows:  . 

1.  Interference  effects  which  involve  changes  in  boundary  conditions,  and  lead  to  iter¬ 
ations  are  not  considered. 

2.  The  effects  of  various  singularities  can  be  combined.  Exact  linear  superposition  of 
flows  is,  of  course,  limited  to  those  flows  which  satisfy  the  Laplace  equation.  It  is  assumed 
that  the  effect  of  the  propeller  on  the  foil  is  mainly  potential  in  origin,  whereas  the  effect  of 
the  foil  on  the  propeller  is  essentially  viscous  in  origin.  The  effect  on  the  propeller  is  treated 
however,  as  a  potential  problem  through  the  use  of  the  total  wake  at  the  propeller. 


3.  A  distribution  of  viscous  wake  behind  a  hydrofoil  is  obtained  by  assuming  that  the 
static  pressure  is  zero  within  the  wake.  This  assumption  seems  reasonable  since  the  pres¬ 
sure  decays  rapidly  downstream,  behind  a  well-streamlined  foil. 

4.  It  is  assumed  that  a. propeller  can  be  simulated  by  a  distributed  sink  disk  of  circumfer¬ 
entially  variable  strength.  This  is  an  adequate  representation  for  calculating  thrust  deduction 
(for  stern  propellers)  since  physically  the  auction  Held  ahead  of  a  propeller,  is  important  and 
the  velocity  field  ahead  of  a  propeller  has  no  vorticity. 

5. .  The  assumption  is  made  that  quasi-steady  propeller  forces  can  be  used  to  estimate 
propeller-thrust  fluctuations  due  to  a  nonuniform  wake  distribution. 

'.ANALYSIS-  "  '• . 

EXPRESSIONS  FOR  THE  FOIL-RESISTANCE  AUGMENTATION 

A  body  in  a  fluid  may  experience  a  lift  force  normal  to  the  flow  direction.  The  force 
assuming  potential  flow  is  given  by  the  well-known  Kutta-Joukowsky  relation 

•'  fl  -Pvr 


where  FL  is  the  lift  farce  per  unit  length, 
p  is  the  mass  density  of  the  fluid, 

V  is  the  undisturbed  velocity,  and 
T  is  the  circulation  around  the  body. . 

Vortex  distributions  can  be  used  in  combination  with  other  flows  to  represent  solid  bodies. 
Flows  about  hydrofoils  with  circulation  can  be  generated  in  this  way. 

Lagally’s  steady-motion  equation5 . 


connects  the  force  Fg  experienced  by  a  source  at  some  point  in  the  flow  field  with  the  output  . 

Qi  '  ••  •  ■. 

Ql  and  the  flow  velocity - at  the  location  of  the  source  QXt  where  r  is  the  distance  from 

4W2  '  \  !  •• 

a  source  Q7  to  a  source  Qx.  When  the  single  sources  of  Equation  [1]  are  replaced  by  line- 
and  surface -source  distributions,  the  interaction  force  (drag  augment)  of  a  propeller-hydrofoil 
system  is  obtained  by  using  the  disturbing  velocities 


WPV  = 


/ 


q(x)dx 

irrr2 
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from  the  foil  acting  at  the  surface-sink  distribution  (propeller)  of  input 


JC 


<7*  (A)  dA 


Wp  is  the  absolute  nondimens ional  perturbation  velocity  which  represents  the  disturbance  of 
the  velocity  V  caused  by  the  foil  in  a  potential  flow.  This,  approach  is  a  matter  of  convenience 
inasmuch  as  the  same  force  would  arise  by  considering  the  propeller-induced  velocity  field  as 
acting  at  the  foil.  To  be  equivalent,  however,  it  should  be  noted  that  only  the  potential  part 
of  the  foil  velocity  is  used  at  the  propeller.  The  arrangement  and  coordinates  of  the  propeller- 


coefficient  CD  and  die  mean-drag* coefficient  CD  Equation  [1]  can  be  written  in  terms  of  the 
foil-distprbing potential  velocities  at  the  propeller  disk  and  a  distributed  surface  sink  of 
total  strength  . 

Dimensionally,  the  relation  is 


D  -  pV  f  <J*WprdrdO 
•'O  •'hub 


[2] 
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and  the  total  augmented  drag  is  given  as  a  coefficient  by  writing  Equation  [2]  in  the  following 
nondimensional  form: 


Cn  =  ■ 


D 


PV2(1-W)R2 


-~r  WpxdxdO 


t3] 


The  total  surface-sink  input  q*  is  related  to  propeller  total  thrust-loading  coefficient  from 
momentum  considerations  by 

?*  =  [-l  +  (l  +  CTf]  Va  [4] 

Since  the  propeller  thrust  varies  with  propeller-blade  position,  the  sink  input  q*  also  varies 
circumferentially.  Similar  to  Equation  [4],  a  mean  surface-sink  input  qf  is  related  to  propeller 
mean  thrust  with 

fff-[-I  +  (l  +  Cr)»]Fa  [4a] 

where  C  T  is  the  mean  thrust-loading  coefficient,  and  a  mean  augmented -drag  coefficient  is 
given  by 


D, 


On  = 


pV2(l-W)R 2 


V~a 


WpxdxdO 


[5] 


Interaction  force  has  been  discussed  as  an  augmented  drag  D  or  augmented-drag 
coefficient  CD.  The  expression  D  =  T-R  (where  T  is  propeller- fluctuating  total  thrust  and 
R  is  foil  resistance  without  a  propeller)  is  variously  called  the  thrust-deduction  force  or 
resistance  augmentation.  The  ratio  of  the  difference  D  to  propeller  average  thrust  T  is 
known  as  the  thrust-deduction  coefficient  t.  The  thrust-deduction  coefficient  can  be  related 
to  the  augmented-drag  coefficient  by 


t 


D  2CD 

T  CTn(l-W) 


T  CTn(l-W) 


EXPRESSIONS  FOR  THE  VISCOUS  WAKE  OF  HYDROFOILS 


The  National  Advisory  Committee  for  Aeronautics  has  made  theoretical  studies  and 
detailed  measurements  of  the  wake  behind  various  airfoils  and  has  derived  empirical  equa¬ 
tions  from  some  of  these  data.8,9  Moreover,  present  knowledge  concerning  free  turbulent 
shear  flows  has  established  the  relevant  wake  characteristics,  e.g.,  the  total  pressure  loss 


5 


is  a  maximum  in  the  center  of  the  wake  and  decreases  to  zero  at  the  wake  edge.  Turbulent 
mixingcauses  the  total  pressure  loss  to  decrease  with  distance  downstream  and  a  consequent 
increase  in  wake  width.  Some  pertinent  empirical  equations,  which  resulted  from  the  subject 
method  of  analysis,  follow. 

Wake  width,  and  maximum  loss  of  dynamic  pressure  in  the  wake  is  given  in  terms  of 
profile-drag  coefficient  and  distance  behind  the  trailing  edge  of  the  foil.  The  wake  half-width 
is 

c  -  0.88  (£+0.15)*  [6] 


where  £  is  the  wake  half-width  in  chords, 

£  is  the  distance  behind  the  foil  trailing  edge  in  chords,  and 
Cio  is  the  section  profile-drag  coefficient. 


The  maximum  loss  of  dynamic  pressure  .is 

2.42  C?0 
£+0.3  * 


[71 


where  rf  is  the  loss  in  dynamic  pressure  at  the  wake  center  as  a  fraction  of  q0.  A  distribution 
of  dynamic  pressure  within  the  wake  is  given  as.  a  cosine  squared  function 


It  is  assumed  that,  within  the  wake,  the  Static  pressure  is  zero  and,  at  the  edge  of  the  wake, 
the  velocity  is  everwhere,  freesteam  velocity  V.  Based  on  the  foregoing,  the  wake  factor  is 
estimated  from  the  relation  . 


<1 


1  - 


2.42 

£+0.8 


It  is  seen  in  Equation  [7]  that  maximum  loss  of  dynamic  pressure  varies  inversely  as  the  first 
power  of  distance  £  downstream.  For  our  purpose,  the  most  important  limitation  of  estimating 
the  wake  factor  from  Equation  [9]  is  the  fact  that  the  relation  was  derived  from  experiments 
with' finite  wings.  However,  aspect  ratios  for  the  test  wings*  range  from  about  3  to  10. 

If  the  wake  centerline  is  assumed  to  coincide  with  the  centerline  of  the  trailing  vortex  • 
sheet,  the  wake  centerline  at  any  given  longitudinal  position  can  be  located  by  the  empirical 
equation10  .  . 


d 

c 


CL  -0.2^/+  0.025 


[10] 
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Relative  to  Ilydrofoil 

The  arrangement  and  quantities  are  shown  schematically  in  Figure  2.  The  empirical  constant  f 
which  is  tabulated  in  Figure  2  was  derived  in  Reference  10  from  studies  of  numerous'  wakes. 

PROPELLER-THRUST  FLUCTUATIONS 

It  is  assumed  that  propeller-thrust  fluctuations  due  to  a  nonuniform  wake  distribution 
over  the  propeller  disk  can  be  estimated  by  calculating  quasi-steady  propeller  forces.  An 
adaptation  of  Burrill’s  strip-rtheory  method  of  calculation  for  marine  propellers  is  used  here 
to  calculate  these  quasi-steady  state  forces.11,12  In  this  method,  unsteady  effects  are  not 
considered;  , i.e.,  when  the  frequency  of  oscillating  velocities  is  small,  a  rotating  propeller 
in  a  variable  wake  is  assumed  to  experience  lift  and  drag  at  each  blade  element  as  if  the 
blade  were  in  steady  flow  at  each  instant.  The  various  angles  and  velocities  at  a  typical 
propeller-blade  section  are  given  in  Figure  3.  When  placed  in  the  wake  of  a  hydrofoil  with 
lift,  a  propeller  will  have  a  down  wash  component  added  to  the  section  inflow  velocity.  In 
general,  variations  per  revolution,  in  propeller  total  thrust  are  essentially  zero  due  to  uniform 
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downwash  because  only  nonlinear  effects  acting  on  the  section  characteristics  contribute  to 
the  thrust  variation  per  revolution.  For  the  most  part,  section  characteristics  may  be  obtained 
from  Reference  11. 


WAKE 


COMPUTED  EXAMPLE 


A  nominal  value  of  the  drag  coefficient  (C^0  =  0.01)*  was  used  to  compute  a  wake  distri¬ 
bution  at  several  positions  behind  a  NACA  16-309  foil.13  The  augmented  resistance  of  the  sub¬ 
ject  foil  was  obtained  for  C ^  values  of  0  and  0.3  (design  lift).  In  this  range  of  lift  coefficient 
(al0  m  -2.04  deg  to  a®  0.92  deg),  the  effect  of  section-drag  coefficient  on  the  wake  distribution 
is  not  too  important  because  of  the  small  changes  in  Cj0  that  occur  in  this  region.13  Likewise, 
the  downward  displacement  of  the  wake  centerline  at  design  lift  is  calculated  from  Equation  [10] 
to  be  negligible  at  £  »  0.15  and  0.30.  Wake  profiles  as  computed  from  Equations  [6]  and  [9]  are 
depicted  in  Figure  4  for  several  positions  downstream  from  the  16-309  hydrofoil.  These  results 
are  used  in  the  calculation  of  propeller-fluctuating  thrust. 


Figure  4  -  Wake  Profiles  at  Several  Positions  behind 
a  NACA  16-3  09  Hydrofoil 


•The  nominal  value  <101  is  confirmed  by  the  resistance-test  results  in  Figure  12.  In  Figure  12,  the  wetted 
surface  is  twice  ths  plan-form  area  and,  therefore,  CjQ  =  2 Ct. 


Disturbing  potential  velocities  Wp  V  from  the  16-309  hydrofoil  were  obtained  from  a 
high-speed  computer  program  of  the  Douglas  Aircraft  Company,  El  Segundo,  California,  and 
is  known  at  the  Taylor  Model  Basin  as  the  Douglas  program.  In  the  method  of  calculating 
potential  flows  with  circulation  as  developed  by  Douglas,  a  solution  is  secured  by  specify¬ 
ing  a  combination  of  a  rectilinear  flow,  a  crossflow,  and  a  vortex  flow.  The  two-dimensional 
potential  flow  behind  the  l«-309  hydrofoil  is  presented  graphically  in  Appendix  A  for  CL  =  0 
and  >=0.3  at  several  distances  (£)  downstream.  The  ordinates  of  the  curves  in  these 
figures  are  the  potential  wake  fraction  Wp  at  the  propeller.  Following  D.W.  Taylor,  a  general 

y  -  Va 

wake  fraction  is  defined  as  W  =  —  . 


PROPELLER-FLUCTUATING  THRUST 


Computational  results  for  fluctuating  total  thrust  are  presented  graphically  in  Figure  5 
for  the  following  propeller  positions  behind  the  NACA  16-309  hydrofoil:  Foil  on  propeller- 


Proptilar  Blotfa  Petition,  In  fegran 


Figure  5  —  Curves  of  Fluctuating  Total  Thrust  for  Propeller  3834 


shaft  axis  with  £  =  0,  0.15,  and  0.30;  and  foil  off  the  propeller-shaft  axis  at  x  =  0.7  for  0.15. 
It  is  noted  that,  when  the  shaft  axis  of  a  3-bladed  propeller  is  parallel  to  the  direction  of  motion 
and  coincides  with  a  symmetrical  wake  centerline,  the  period  is  60  deg  because  of  symmetry. 

360 

If  a  3-bladed  propeller  is  off  center,  the  period  is  - =  120  deg.  These  total  thrust  curves 

3 

are  used  subsequently  to  calculate  interaction  forces. 

A  harmonic  analysis  of  the  fluctuating  thrust  curves  was  made  using  a  12-ordinate 
scheme  given  by  Scarborough.14  The  dimensional  Euler  coefficients  an  and  bn  are  given  in 
Table  1  for  each  propeller-hydrofoil  arrangement.  It  is  noted  that  only  cosine  terms  appear  in 


TABLE  1 


Harmonic  Coefficients  in  the  Series  Expansion 


Foil  on  Propeller  Axis* 

Fori  off  Axis  at  ,r  =  0.7**  . 

Harmonic  Coefficients 

4T=  o.oo 

£=0.15 

•£=0.30 

. 

CP 

'll 

% 

46.2 

46.6 

47.3 

53.37 

10.66 

10.28 

9.57 

0 

’  ^  2 

CO 

b> 

cr» 

:  6.73 

5. -5  8 

1.620 

a3 

6.52 

4.03 

2.58  ■ 

»  o 

4.98 

2.50  . 

0.683 

-  0.487 

4.17 

1.68 

-  0.303  ' 

.o' 

1.97 

0.666 

-  0.308  ' 

■  -  0.589 

h 

0 

0 

n 

\j 

-  5.17 

'  h2 

0 

0 

0 

.0 

0 

■  0 

■  0 

-  0.175 

0 

0 

0  . 

0 

h  . 

0 

0 

0 

0:507 

•Period  =  60  degrees 
♦•Period  =  120  degrees 


the  expansion  for  the  foil  on  the  propeller  shaft  axis  for  £^0.3.  This  is  caused  by  the  symmetry 
of  the  cosine-squared  wake  distribution  and  its  location.  For  the  off-axis  position,  sine  and 
cosine  terms  alternate.  The  Euler  coefficients  can  be  combined  as 

A„  sin  («?♦  <?„) 


where 


4  =  \A*2  +  bn 


<f>n  =  tan' 


■a 


Plotted  in  Figure  6  as  a  function  of  £  are  the  resultant  positive  harmonic  coefficients  An  expressed 
as  a  fraction  of  the  direct  component  for  each  on-axis  condition.  When  the  propeller  has  reached 
the  £  *  0.30  position  downstream,  apparently  only  the  first,  second,  and  third  resultant  harmonics 
contribute  significantly  to  total  thrust  fluctuation.  The  coefficients  An  for  the  foil  off  the 
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Figure  6  -  Variation  of  Resultant  i'arrnonic  Coefficients  of 
Fluctuating  Thrust  with  £ 

Curves  arc  for  foil  on  propeller-shaft  axis. 


propeller  axis  ( x  =  0.7)  are: 


£  =  0.15 


Order 

c0 

1 

0.0270 

?.  •  ! 

0.0304 

3  ; 

0.0033 

i 

0.0091 

5 

0.0095 

6 

0.0110 

THRUST  DEDUCTION 

Equations  [3]  and  [5]  are  evaluated  numerically  using  Simpson’s  rule.  The  functions 
<7*/Fa  and  7f/Fa  are  calculated  using  the  appropriate  total-thrust  curve  in  Figure  5.  In  cal¬ 
culating  Vi/Va,  the  direct  component  may  be  obtained  directly  from  Table  1.  Appropriate 
values  of  Wp  are  read  from  the  curves  in  Appendix  A  and,  with  the  input  data  as  outlined, 
numerical  solutions  were  obtained  for  several  assumed  values  of  effective  wake  factor  (1  -  W) 
Resulting  values  of  and  =  Cfi  -  Cj)^  are  tabulated  in  Table  2  for  several 

propeller-hydrofoil  arrangements  and  conditions. 
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TABLE  2 


Computed  Augmented -Drag  Coefficients  for  Various 
Propeller-Hydrofoil  Arrangements  and  Conditions 


e 

Foil 

Lift 

cL  =  o.o 

CL  =  0.3 

(1-1V) 

10  2  cn 

io  2cDl 

102  CD 
u2 

102Ct) 

io2c0i 

102 

V 

0 

EH 

nm 

1.645 

1.731 

m 

1.284 

1.386 

0.97 

- 

- 

- 

1.191 

1.074 

0.15 

0.80 

1.026 

1.044 

0 

1.043 

11P1B 

0.90 

0.8204 

0.8340 

!=;  0 

0.8334 

iS|s« 

0.97 

- 

- 

- 

0.7151 

0.30 

0.80 

■■pppi 

0.6833 

0 

0.6701 

0.6658 

0.90 

ES 

0.5458 

~0 

0.5353 

0.5373 

0.97 

- 

- 

- 

0.4581 

0.4652 

Foil  off  Axis  at 

r  =  0.7 

0.15 

1.340 

1.348 

~  0 

1.071 

1.078 

~  0 

0.97 

0.9204 

0.9253 

—  0 

In  Figures  7  and  8,  the  thrust-deduction  curves  were  computed  using  the  augmentation 
coefficients  of  Table  2  and  the  relations 

2C0  2(?01 

t  - -  ,  - - 

C'Tn(l-W)  CTn(l-W) 

The  effects  of  propeller-hydrofoil  spacing,  wake,  fluctuating  thrust,  and  foil  lift  are  all  ap¬ 
parent  in  these  figures.  For  (1  -  IV)  =  0.90,  which  is  more  likely  to  be  obtained  in  practice, 
the  effect  on  t  and  tx  of  going  from  =  0  to  =  0.3  is  not  important  because  the  effect  on 
total  system  (1  -  t)  is  quite  small  even  in  the  vicinity  of  £=  0. 

Considering  conditions  at  design-lift  coefficient  (C L  =  0.3),  the  influence  of  spacing  f 
on  the  thrust-deduction  coefficients  t  and  is  significant.  In  going  from  £  =  0  to  £  =  0.3,  t 
and  tx  are  reduced  from  an  order  of  magnitude  of  approximately  0.05  to  approximately  0.02. 

The  higher  the  effective  wake  W ,  the  more  pronounced  is  this  reduction.  The  contribution  of 
the  fluctuating  part  of  the  total  propeller  thrust  on  t  is  apparent  only  at  or  near  the  foil  trailing 
edge  (£=  0).  Relative  to  the  on-axis  condition,  the  off-axis  curve  in  Figure  8  appears  parallel 
and  displaced  upwards  about  0.005  units.  To  further  investigate  off-axis  performance,  addi¬ 
tional  calculations  were  made  at  x  =  0.4  (in  lift  direction)  based  on  an  average  thrust-loading 
coefficient  C T  =  0.1979.  The  average  thrust  coefficient  is  based  on  the  direct  thrust  compo¬ 
nents  given  in  Table  1,  and  an  assumed  value  of  0.90  for  the  wake  factor  (validated  later  by 
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Figure  7  —  Computed  Thrust-Deduction  Coefficient  versus  £for 
Hydrofoil  16-800  and  Propeller  3834 


Figure  8  -  Computed  Thrust-Deduction 
Coefficient  versus  (l-W')  at  C ^  =  0.3 
for  Hydrofoil  16-309  and 
Propeller  3834 


the  experimental  results  in  Table  4).  A  computed  value  of  t  =  0.036  was  obtained  for  x  =  0.4. 
This  result  is  the  same  as  that  obtained  for  the  off-axis  condition  at  x  =  0,7  (see  Figure  8). 

The  effective  wake  is  essentially  constant  within  the  range  of  *  =  0.4  to  0.7;  consequently, 
the  average  thrust  is  the  same  at  the  same  revolutions.  This  means  that  the  same  sink  strength 
combines  with  a  weaker  disturbing  field  Wp  to  give  a  lower  thrust  deduction.  The  accuracy  to 
which  this  weaker  potential  field  can  be  determined  is  questionable.  In  any  event,  a  theoreti¬ 
cally  obtained  curve  of  thrust  deduction  with  off-axis  position  would  tend  to  zero  for  foil  loca¬ 
tions  out  of  the  propeller  disk. 


EXPERIMENTS 

DESCRIPTION  OF  PROPELLER  AND  HYDROFOIL 

The  TMB  stock  propeller  selected  for  the  propeller-hydrofoil  configuration  was  3834, 
a  three-bladed  aluminum  propeller  with  a  diameter  of  8  in.  A  drawing  of  this  propeller  is 
shown  in  Figure  9.  The  NACA  16-309  high-lift  and  low -drag  section  was  used  for  the  hydro¬ 
foil  (TMB  Model  H-41)  which  corresponds  to  the  submerged  hydrofoils  of  Navy’s  first  opera¬ 
tional  hydrofoil  patrol  boat  PC(H).  The  offsets  for  the  NACA  16-309  section  are  tabulated 
in  Table  3.  The  ratios  foil  chord  to  propeller  diameter  c/d  and  foil  maximum  thickness  to 
propeller  diameter  t/d  are  1.67  and  0.15,  respectively. 

Model  H-41  was  equipped  with  a  removable  nacelle  of  3.5-in.  diameter  (TMB  Series-58 
shape)  and  an  adjustable  tail  flap.  Provision  was  made  for  removing  segments  of  the  foil  from 


the  lower  end.  With  these  features  incorporated  in  the  hydrofoil,  the  augmented  resistance 
could  be  determined  for  test  conditions  which  represent,  with  and  without  nacelle,  infinite 
span  and  semi-infinite  span.  A  parameter  designated  the  span  ratio  l/L  was  used  to  define 
the  geometry  of  the  semi-infinite  span  conditions.  In  the  span  ratio,  l  is  the  foil  span  meas¬ 
ured  from  the  axis  of  the  water  tunnel  jet  to  the  free  end  of  the  foil  and  L  is  the  foil  span 
measured  from  the  top  of  the  water  tunnel  jet  to  the  free  end  of  the  foil. 

ARRANGEMENT  AND  INSTRUMENTATION  FOR  TESTS 

Figure  10  shows  hydrofoil  Model  11-41  and  the  drag  dynamometer  as  assembled  prior  to 
conducting  tests  in  the  24-in.  water  tunnel.  Figure  11  shows  the  arrangement  for  a  typical 
test  of  the  hydrofoil  and  propeller  in  this  tunnel.  The  foil  was  mounted  vertically  across  the 
tunnel  jet.  A  column  which  was  secured  to  a  drag  dynamometer  was  connected  (outside  the 
tunnel  jet)  to  the  upper  end  of  the  foil.  The  lower  end  of  the  foil  extended  beyond  the  jet 
boundary  and  was  free.  The  hydrofoil  and  propeller  were  not  mechanically  connected  so  the 
foil-resistance  augment  is  obtained  directly  from  the  resistance  measurements. 

The  drag  dynamometer  was  of  the  TMB  magnigage  type  and  had  a  linear  calibration 
w'ithin  the  range  of  drag  and  lift  forces  encountered  in  the  investigation.  The  accuracy  of  the 
drag  measurements  was  to  ±1  percent  of  full  load.  Accuracy  of  propeller  measurements, 
using  the  10-hp  dynamometer  for  the  24-in.  water  tunnel  was  also  about  ±1  percent  of  full 
load.  • 

TEST  PROCEDURE 

Resistance  tests  with  and  without  the  propeller  were  conducted  at  atmospheric  pressure 
for  each  hydrofoil  configuration  and  position.  All  tests  were  conducted  at  a  constant  speed  of 
14  knots  except  Test  4  (flap-angle  variation)  which  was  run  at  a  speed  of  12  knots  because  ex¬ 
cessive  vibration  was  encountered  at  14  knots.  The  constant-speed  tests  with  the  propeller 
were  conducted  by  measuring  propeller  thrust  and  torque  and  foil  resistance  over  a  range  of 
propeller  rpm  which  corresponded  to  J  values  of  about  0.7  to  1.1.  In  addition  to  constant- 
speed  runs,  the  effect  of  test  Reynolds  number  (speeds  of  6  to  14  knots)  on  the  hydrofoil  re¬ 
sistance  for  Test  1  was  investigated  for  the  purpose  of  establishing  the  adequacy  of  test 
Reynolds  number.  ^ 

DISCUSSION  OF  RESULTS 

Thrust-deduction  coefficients  versus  propeller-speed  coefficient  are  given  in  Appen¬ 
dix  B  for  each  test.  Experimental  data  points  are  plotted  for  each  curve.  Although  the  results 
are  discussed  primarily  in  terms  of  faired  thrust-deduction  coefficients  for  a  propeller -speed 
coefficient  J  =  0.9  (corresponds  closely  to  maximum  propeller  efficiency),  a  wide  range  of  J’s 
was  covered  in  order  to  provide  additional  data  for  fairing  purposes  and  to  investigate  the 
effect  of  J  on  t  more  completely. 
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Figure  11a  —  Test  3  with  Foil  at  0.30  Position 


Figure  10  -  TMB  Hydrofoil  Model  H-41  and 
.  Dynamometer  Assembly 


Figure  lib  —  Test  3  with  Foil  at  0.15  Position 

Figure  11  -  Model  H-41  Installed  for  Testing 
in  24-Inch  Water  Tunnel 
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The  results  of  tests  to  determine  Reynolds  effect  are  presented  in  Figure  12  as  curves 
of  total -resistance  coefficient  Ct  versus  Reynolds  number  Rn  with 


and 


Rt 


'ApSV2 


It  is  seen  that,  for  £  =  0.15  and  £=  0.30,  the  Ct  curves  appear  to  parallel  the  Schoenherr 
turbulent  frictional-resistance  line  for  a  flat  plate. 

Table  4  summarizes  the  test  conditions  and  corresponding  experimental  values  of 
thrust-deduction  coefficient  t  and  wake  factor  (l-lV).  A  comparison  of  the  experimental 
values  of  (i-fT)  for  Test  2  and  the  computed  on-axis  data  which  are  cross-plotted  in 
Figure  8  shows  that  the  t  values  agreed  for  both  the  £=  0.15  and  £=  0.30  positions: 


Thrust-Deduction  Coefficient 

£ 

(l-W) 

Theoretical 

Experimental 

0.15 

0.910 

0.031 

0.031 

0.30 

0.919 

0.019 

0.018 

Examination  of  the  off-axis  experimental  results  and  the  theoretical  results  that  are 
plotted  in  Figure  13  indicates  that  the  experimental  curve  showed  a  slight  decrease  in  t 
(from  0.031  to  0.021  in  going  from  x  =  0  to  x  =  0.7  in  the  lift  direction).  This  reduction  in 
t  is  to  be  expected  and,  of  course,  t  becomes  zero  when  the  foil  is  moved  completely  away 
from  the  propeller  field.  Although  the  theoretically  obtained  thrust  deduction  remained  es¬ 
sentially  constant  at  t  =  0.036  for  a  range  of  x  comparable  to  the  experiments,  the  order  of 
magnitude  for  t  was  correct. 

The  most  important  result  shown  by  the  remaining  test  data  in  Table  4  and  in  Figures  14 
and  15  is  the  small  change  in  one-minus-the-thrust  deduction  (a  negligible  effect  on  the  total 
system  due  to  thrust-deduction  variations)  for  the  various  configurations.  This  result  confirms, 
of  course,  the  analytical  procedure.  Figure  15  shows  a  negligible  change  in  t  from  l/L  =  0.3 
to  l/L  =  0.5  which  approximates  an  infinite  aspect  ratio.  Since  the  span  ratio  was  varied  by 
removing  segments  of  the  span  from  the  free  end  only,  the  increase  in  foil  drag  due  to  the  flow 
at  one  tip  was  about  one-half  of  that  which  would  occur  for  a  usual  finite  wing. 

The  following  additional  remarks  are  made  relative  to  the  experimental  data: 

1.  The  experimental  results  confirmed  the  theoretical  result  which  indicated  that  the 
effect  of  foil  lift  (CL  maximum  =  0.3)  on  thrust  deduction  can  be  disregarded  except  at  or 
near  zero  clearances. 
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Thrust  Deduction  Coefficient 


Figure  12  —  Total-Resistance  Coefficient  versus  Reynolds 
Number  for  Model  B*ll,  Test  1 
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Figure  13  -  Variation  of  Thrust-Deduction  Coefficient  with  Off-Axis  Position, 

Tests  2  and  6 
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TABLE  4 


^  (0.7/?  in  lift  direction) 

0.86  0.3  0.15  0.50  0  X  0.032  0.913  Off  Axis 

(0.66/?  opposite  to  lift  direction) 


Thrust- Deduction  Coefficient,  t 


Model  H-41 
Propeller  3834 
Test  3  and  4 
♦  *  0.15 
1/ L*  0.50 


\)  =  0.9 


I  23456789  10  II 

Plop  Angle  in  degrees 


Figure  14  -  Thrust-Deduction  Coefficient  versus  Flap  Angle  for  Foil  with  Nacelle 

Tests  3  and  4 


Model  H-41 
Propeller  3034 
Test  3  and  5 
a*-0. 20° 

«  =0.15 


t  for  J*0.90 


Spon  Rotio  /(_ 

Figure  15  -  Variation  of  Thrust-Deduction  Coefficient  and  Foil  Drag  per  Unit 
Length  with  Span  Ratio  for  Foil  with  Nacelle,  Tests  3  and  5 


2.  Within  the  normal  operating  range  of  the  propeller,  the  effect  of  propeller-speed  coef¬ 
ficient  on  the  thrust-deduction  coefficient  was  nil. 

CONCLUDING  REMARKS 

A  method  which  is  based  on  Lagally’s  steady-motion  theorem  was  derived  and  used  to 
calculate  the  resistance  augmentation  for  a  submerged  hydrofoil.  In  the  method,  mathematical 
singularities  are  derived  from  propeller-fluctuating  thrusts  as  estimated  from  nuasi-steadv 
theory  using  empirically  obtained  wake  data.  The  disturbing  foil  velocities  for  potential  flow 
are  obtained  from  a  high-speed  computer  program.  A  harmonic  analysis  of  the  total  fluctuating 
thrust  curves  computed  for  a  three-bladed  propeller  indicated  that,  when  the  propeller  was 
positioned  about  0.3  times  the  foil  chord  downstream  from  the  trailing  edge  and  on  an  axis 
through  the  foil  trailing  edge,  only  the  first,  second,  and  third  resultant  harmonics  contributed 
significantly  to  total  thrust  fluctuation.  Computed  and  experimental  thrust-deduction  results 
for  several  hypothetical  on-axis  and  off-axis  arrangements  and  conditions  were  in  good 
agreement. 

The  most  salient  features  concerning  the  effects  on  thrust  deduction  of  the  parameters 
and  conditions  investigated  are  summarized  as  follows: 

1.  The  contribution  of  the  fluctuating  part  of  the  propeller  total  thrust  to  the  thrust  deduc¬ 
tion  was  important  when  the  propeller  was  located  at  or  near  the  foil  trailing  edge. 

2.  Clearance  between  the  foil  trailing  edge  and  the  propeller  was  important.  In  going  from 
zero  clearance  to  0.3-chord  length,  one-minus-the-thrust  deduction  varied  from  0.916  to  0.980. 

3.  In  the  range  of  lift  coefficients  from  zero  to  design  lift  ( C L  =  0.3),  the  influence  of  lift 
on  one-minus-the-thrust  deduction  was  generally  negligible. 

4.  Within  the  normal  operating  range  of  the  propeller,  the  effect  of  propeller-speed  coef¬ 
ficient  on  thrust  deduction  was  nil. 

5.  Within  the  range  of  investigation,  there  was  little  effect  on  the  total  system  due  to  thrust- 
deduction  variations  with  foil-span  ratio  (0.2  8  to  0.50)  and  foil  flap  angle  (0  to  10  deg). 
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Figure  16  —  Computed  Potential  Wake  Fraction,  Foil  16-309,  CL  =  0,  <*=-2.05  Degrees 
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Figure  18c  -  Foil  Off  Axis  at  0.9R  in  Lift  Direction 
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Propeller  Speed  Coafficiant  J 

Figure  19  -  Thrust-Deduction  Coefficient  versus  Propeller-Speed  Coefficient  with  £  and  CL  as 
Parameters  for  Foil  without  Nacelle,  Tests  1  and  2 
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Figure  20  -  Thrust-Deduction  Coefficient  versus  Propeller-Speed  Coefficient,  with  f  and  CL  as 

Parameters  for  Foil  with  Nacelle,  Test  3, 
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Figure  21  —  Thrust-Deduction  Coefficient  versus  Propeller-Speed  Coefficient  with  Flap  Angle  as  a 

Parameter  for  Foil  with  Nacelle,  Test  4 
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Figure  22  -  Thrust-Deduction  Coefficient  versus  Propeller-Speed  Coefficient  with  l/L  as  a 

'Parameter  for  Foil  with  Nacelle,  Test  5 
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